Maximum relative excitation of a specific vibrational mode 
via optimum laser pulse duration 
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For molecules and materials responding to femtosecond-scale optical laser pulses, we predict max- 
imum relative excitation of a Raman-active vibrational mode with period T when the pulse has an 
FWHM duration r ~ 0.42 T. This result follows from a general analytical model, and is precisely 
confirmed by detailed density-functional-based dynamical simulations for Ceo and a carbon nan- 
otube, which include anharmonicity, nonlinearity, no assumptions about the polarizability tensor, 
and no averaging over rapid oscillations within the pulse. The mode specificity is, of course, best 
at low temperature and for pulses that are electronically off-resonance, and the energy deposited in 
any mode is proportional to the fourth power of the electric field. 



For a quarter century there has been considerable in- 
terest in optimizing the vibrational response of molecules 
and materials to ultrafast laser pulses |ll4l5[. This prob- 
lem is directly relevant to the broader issue of coherent 
control in physical, chemical [l6l - [l9| . and biological (20I - 
[23} systems. 

Here we consider excitation via impulsive stimu- 
lated Raman scattering and related techniqties using 
femtosecond-scale optical pulses. We find that the opti- 
mum fuU-width-at-half-maximum (FWHM) pulse dura- 
tion for exciting a specific vibrational mode with angular 
frequency 2it/T is given by r ~ 0.42 T. Our prediction 
results from a general analytical model, and is precisely 
confirmed by completely independent density-functional- 
based simulations for Cgo [2J-|28| and a small carbon 
nanotube (29l - [3H . Unlike the model, these simulations 
include anharmonic effects in the vibrations, nonlinear 
effects in the response to the applied field, and no sim- 
plifying assumptions about the electronic polarizability 
tensor. 

Our general model consists of the following: 
(1) The electric field has the form 



£ {t) = £q sin (7rt/2T) sin {ujt + S) 



< i < 2r (1) 



with w 3> tt/t and huj off resonance. Since the oscilla- 
tions of sin^ {ujt + 6) will average out to 1/2 over a period 
that is short compared to the response time of the vibrat- 
ing nuclei, we will actually replace the square of Eq. ([T]) 
by the envelope function 



£{tf=Po 



(7rt/2T) 



< t < 2r 



(2) 



with £n_r^ ^o/2- This form has the following nice fea- 
tures [32[: (i) The duration is finite and need not be 
truncated, (ii) The FWHM duration is exactly half the 
full duration 2t. (iii) A plot reveals that it closely resem- 
bles a Gaussian, (iv) The slope is zero at beginning and 
end. 



(2) The initial conditions Qk (0) = Qk (0) = are 
imposed on the normal-mode coordinates Qk (t). This 
approximation is valid for the expectation value below 
Eq. (|4]) at low temperature, or an ensemble average at 
higher temperature, in the linearized Eq. ([S]). 

(3) The equation of motion is given by the standard 
(Placzek) model for Raman-active modes [1, [2I, 33 1: 



d^Qk/dt^ 



{t f /2 



(3) 



where the anharmonic and damping terms (in the 
normal-mode coordinates Qk), nonlinear terms (in the 
electric field £ (t)), and off-diagonal terms (in the polar- 
izability tensor a) have been neglected, with the electric 
dipole moment given by /xtot = M+Mind, Mind = «•£+..., 
<^xx = ao + J2k '^'kQk + ■ • and the radiation field taken 
to be linearly polarized in the cc-direction. Equation (|3]) 
follows from e.g. the Heisenberg equations of motion for 
the quantized Hamiltonian 

^ = E (^fe /2 + ^iQl/^) -(^■£-£ S,-£/2 (4) 

A; 

with Qk = (Qk'^ and aj^ = (^da/dQkJ, when the term 

involving dfijdQk vanishes or is neglected (e.g. in taking 
rapid oscillations of the electric field £ (t) to average out 
to zero, an approximation which also causes the lowest- 
order nonlinear term to vanish in Eq. (|3])). 

If Eq. ([2]) is substituted into Eq. (|3]), the solution after 
the pulse is found to reduce to 



Qk{t) 

where 



a^f^ sin(rfe7r) 



4-^ 



rk 



sin(cjfet — r^Tr) , t > 2t (5) 



2r/Tfc 



(6) 
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The total energy (^Q^, + i^f-Q^ /2 in vibrational mode k 
is therefore 



El, 



1 

32 



a'f, sin (rfcTr) 



c4 
^0 



t > 2t 



(7) 



i?fe is, of course, equivalent to the maximum kinetic en- 
ergy if™'"''. 

The maximum response is then given by the extremum 
of the function in brackets, which occurs at rk ~ 0.8375, 
or 2t « 0.8375 Tfe. 

We have tested this prediction by performing inde- 
pendent supercomputer simulations for Cgo and a small 
carbon nanotube, using the density-functional-based 

"""" " to- 



M |3£ 



approach of Frauenheim and co-workers 
gether with semiclassical electron-radiation-ion dynamics 
(SERID), which is defined by the following equations [SGf : 
(1) Time-dependent Schrodinger equation in a 
nonorthogonal basis: 



ihdxp (t) /dt= S'^ ■ H - tpit) 



(8) 



With 60 atoms and a minimal basis set, the matrices are 
240 X 240. A time step of 50 attoseconds was used. The 
simulation time after completion of the pulse was 2000 fs 
for Ceo and 1000 fs for the nanotube. 

(2) Ehrenfest's theorem (in a nonorthogonal basis): 



M- 



1 V- ,t fdH ^dS d 



dX dt 



(9) 



where X is any nuclear coordinate. The Hamiltonian 
matrix H , overlap matrix S', and effective ion-ion repul- 
sion Uj.p,r, were determined by the methods and results of 
Refs. [3J] and [s^ and later work by this group. 

(3) Coupling of electrons to the radiation field through 
the time-dependent Peierls substitution 



(10) 



where q = —e. A is the vector potential, which in the 
present simulations was taken to have the form 

A{t) = Aosin{Trt/2T)cos{Lut) , < i < 2t (11) 

with uj S> 7r/r, so that £{t) « lo Aq sin (irt / 2t) sm{ujt), 
although Eq. ifTTj) was actually used. The polarization 
vector was taken to lie along the x-axis, with the z-axis 
pointing down the axis of a nanotube. 

Within the present DFT-based model, the energy gap 
for electronic excitations is 1.80 eV for Cgo and 1.56 eV 
for a (3, 3) nanotube with a periodicity length of 5 unit 
cells. (This model nanotube, with only 60 atoms, has a 
substantial gap because small wavenumbers kz are not al- 
lowed. An infinitely long (3,3) nanotube would be metal- 
lic [iOl.) The laser pulse photon energy was chosen to be 



0.69 eV for Ceo and 0.80 eV for the nanotube, and is thus 
off-resonance in both cases. 

For Ceo, seven of the ten Raman-active modes were 
appreciably excited: modes ^£,(1), Ag{2), Hg{l), Hg{'i), 
Hg{5), Hg{6) and Hg{8), with periods of 59.7 fs, 20.6 fs, 
125 fs, 38.5 fs, 27.0 fs, 23.0 fs, and 19.0 fs respectively. 
For the (3, 3) nanotube only an Eg mode with a period 
of 182 fs was observed, and this is consistent with the 
experimental difficulty of observing the radial breathing 
mode in the Raman spectrum f3l|. 

As in Ref . |12|] , it is natural to characterize the strength 
of the vibrational response of a specific mode by its max- 
imum kinetic energy K^'^^. In the harmonic approxi- 
mation, after completion of the laser pulse, the veloc- 
ity Qk{t) is proportional to cos (cj^t + (5), so is pro- 
portional to cos^ {ujkt -I- 15) = [1 -I- cos {2ojkt + 26)] /2. A 
numerical Fourier transform of the total kinetic energy 
therefore shows a peak at 2wfe, with a strength propor- 
tional to the response of the normal mode with angular 
frequency ujk- Figure [T] shows our results for the Fourier 
transform of the total kinetic energy for Ceo following a 
12 fs, 0.69 eV, 5.0 V/nm pulse. 
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FIG. 1; Numerical Fourier transform of the total kinetic en- 
ergy of Ceo after being subjected to a laser pulse with a 
FWHM duration r of 12 fs, photon energy hu of 0.69 eV, and 
electric field amplitude £o of 5.0 V/nm. Seven Raman-active 
modes are clearly observed, each with a frequency ujk/2-n 
which is half the frequency exhibited in the figure by the cor- 
responding kinetic energy Kk- 

Figure[2]shows that the detailed simulations agree with 
Eq. ([7]) regarding the dependence of the maximum kinetic 
energy iC™^'^ of mode k on the field amplitude £o (for 
the range of amplitudes considered here) when the pulse 
duration is fixed: K^'^^ ex. £q. 

When £o is fixed, on the other hand, K^^^ is predicted 
to be proportional to the square of the function involving 
Tfe in Eq. ([7]), and is maximal when rk ~ 0.84. Figure [3] 
shows, for example, that the Hg{8) mode in Ceo domi- 
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FIG. 2: Maximum kinetic energy of ^9(1), ^,,(2), Hg{l), 
//g(4), Hg{5), Hg{6) axid Hg{8) modes of Ceo responding to 
0.69 eV laser pulses when the FWHM pulse duration r is 
fixed at 12 fs and the electric field amplitude £^0 is varied. 
Here, and in the subsequent figures, each point represents a 
DFT-based SERID simulation and each curve has the form 
predicted by Eq. 0, with no adjustable parameters except 
the effective polarizability parameter a'l^, which is determined 
by a least-squares fit. 



the analytical model (curves) and the completely inde- 
pendent DFT-based simulations (points) is truly remark- 
able. The only adjustable parameter for each curve is the 
effective polarizability parameter a'f.. 

There is similarly remarkable agreement for the carbon 
nanotube, as can be seen in Fig.m In this case, two pho- 
ton energies hoj were used: 0.8 eV and 2.0 eV, which are 
respectively below and above the excitation gap of 1.56 
eV for this model nanotube, as discussed above. Figure|3] 
shows the results for the maximum kinetic energy K^fP^^^ 
of the one prominent mode observed in this case. (Re- 
call that the polarization vector for the pulse was chosen 
to point across the nanotube, so that only radial modes 
should be directly excited, and this is a severe constraint 
for the small (3,3) nanotube, as is also observed in ex- 
periment [sij.) The panels in Fig. m correspond to: (a) 
huj = 0.8 eV, T = 76 fs; (b) huj = 0.8 eV, £0 = 9.66 V/nm; 
(c) huj = 2.0 eV, T = 76 fs; (d) hu = 2.0 eV, £0 ^ 9.66 
V/nm. It is interesting that the analytical model still 
works in this case even for above-the-gap excitation, but 
we find less satisfactory agreement for Ceo, and certainly 
do not expect the model of Eq. ^ to be generally valid 
when there is a substantial population of electrons in ex- 
cited states at the end of the pulse. In this case other 
effects will ordinarily be of dominant importance 0, [HI . 



nates for 2t < 23.4 fs, and the Hg{l) mode for 2r > 23.4 
fs, in agreement with the analytical model of Eqs. dH-pt , 
and in qualitativ e ag reement with other simulations [l2| 



and experiment [2J-|2 
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FIG. 3: Maximum kinetic energy for the 7 prominent modes 
of Ceo responding to 0.69 eV laser pulses when £0 is fixed at 
5.0 V/nm and the full pulse duration time 2r is varied. The 
agreement between the general model of Eqs. ([21)-© (curves) 
and the density-functional-based simulations of Eqs. (f8))- (|ll|) 
(points) is remarkable. 

In Fig. [31 the agreement between the predictions of 
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FIG. 4: Maximum kinetic energy of Eg mode for (3,3) carbon 
nanotube, excited by a 0.8 eV laser pulse (a) when r is fixed at 
76 fs and So is varied and (b) when £o is fixed at 9.66 V/nm 
and T is varied. Panels (c) and (d) are the corresponding 
results for a pulse with a photon energy of 2.0 eV. 

As mentioned above, the electronic response in the 
DFT-based simulations can be characterized by an effec- 
tive polarizability parameter a'^, if one fits Eq. ^ to the 
results of the simulations. If the model leading to Eq. ([7]) 
is in fact consistent with the detailed simulations, then 
the values obtained when (a) varying the electric field 
amplitude £q at constant (FWHM) pulse duration r and 
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TABLE I: For each of the prominent modes observed in the 
response of Ceo and the (3,3) carbon nanotube, the values of 
the effective a'^. were obtained through the procedures desig- 
nated (a) and (b) in the text. In procedure (a), So was varied 
while the (FWHM) duration r was fixed at 12 fs, for Ceo, 
or 76 fs, for the nanotube. In procedure (b), r was varied 
while £q was fixed at 5.0 V/nm, for Ceo, or 9.66 V/nm, for 
the nanotube. The nanotube was also subjected to 2.0 eV 
pulses, and the effective values of aj. obtained through both 
procedures are again shown. 
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(b) varying r at constant £q should also be reasonably 
consistent. Table U shows the results for all of the promi- 
nent Raman-active modes. The two procedures do lead 
to generally consistent values, with the differences for the 
Ag{2)^ Hg[6), and Hg{8) modes presumably arising from 
their short periods, so that the approximation of averag- 
ing out the Lo oscillations in Eq. © is not valid. However, 
the prediction regarding optimum pulse duration works 
even for these modes, as can be seen in Fig. [3l 

It has long been suspected that there might be an op- 
timum pulse duration. 



For example, in Ref. [12| the 
pulse envelope was assumed to have a Gaussian form: 
|f (t)l^ = Aexp (— 2t^/r^) cos^ (uit) (with a slight change 
in notation). Simulations were then performed with a 
model that essentially lies between the analytical model 
of Eqs. (12)- (IS]) and our much more extensive density- 
functional-based simulations, defined by Eqs. (|8|)-([TT|). 
It was found that the optimum value of f was about 
r/3.4. The corresponding FWHM duration then has an 
optimum value r = ^/2\n2 f sa 0.35 T. In both our 
general model and DFT-based simulations, on the other 
hand, the optimal FWHM pulse duration is found to be 
given by T 0.42 T. Both results seem to be consis- 



tent with the existing experiments [2J-|28|, and both are 
larger than the naive prediction of r = 0.25 T. In our 
result, the reason appears to be the phase lag apparent 
in Eq. dS]). 

It should be emphasized that our result of r « 0.42 T 
is for the maximum response of this mode when the in- 
tensity of the laser pulse is kept fixed while the duration 
is varied. On the other hand, if the fluence is instead held 



fixed (while r is varied), Eqs. ([5]) and (O are changed to 

^ , , a'i.W sin(rfc7r) . , , 

Qk (t) = — ^ ' sm Wfet - rkTT) , t>2T 

WfcTTceo rfe(l - r^) 

(12) 

, t>2T (13) 
(14) 



and 



where 



Ek = 



a'^W sin(rfe7r) 



TTceo rfe(l - rl) 



W ~ —-Cnl 



is the energy in the pulse. If W is held constant, the 
maximum in Eq. (jl3p is achieved as rj, — > 0: 



Ek (0) = 



1 



ceo 



1 2 



SO that 



Ekir) 
Ek{0) 



sin (rfeTr) 



7rrfe(l - rl) 



(15) 



(16) 



Then rfe = 0.84 gives Ek (0.42 T) /Ek (0) = 0.385. How- 
ever, Eq. (|15p holds for all modes, so there is no relative 
enhancement of any preferred mode. 

It should also be emphasized that the success of the 
simple analytical model of Eq. ^ results from the fact 
that the additional effects described below Eq. ^ are 
not very important for the pulse intensities, time dura- 
tions, and systems considered here. The various effects 
omitted in the model will lead to richer (or messier) be- 
havior in regimes where nonlinear effects, anharmonic- 
ity, damping, etc. become comparable in importance to 
the lowest-order effects included in the model. In future 
studies with density-functional-based simulations based 
on Eqs. ((8|)- (fTT|l it may be interesting to sort out the 
influence of the higher-order effects. 

Finally, as emphasized below Eq. ([T]), the present re- 
sults are for huj off resonance, in contrast to those of 
Refs. 



e.g. Kets. [IJJJ and [IJJ. In particular. Smith and Cina 
considered pulses whose central frequencies are near res- 
onance with an electronic transition, in a model with 
two electronic levels and one vibrational degree of free- 
dom. They found that the momentum increment (for the 
nuclei) falls off less rapidly with increasing offset (from 
resonance) than does the population loss (from the elec- 
tronic ground state), and they presented a rather com- 
plex strategy for optimizing the various parameters for 
these "preresonant" pulses. The work of this group and 
others is thus complementary to that of the present pa- 
per. 

In summary, we find remarkable agreement between 
the general model of Eqs. and the detailed DFT- 

based simulations based on Eqs. (I5))- ([TT|) for Ceo and a 
carbon nanotube. At fixed pulse intensity, both of these 
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approaches predict maximum excitation of a Raman- 
active vibrational mode with period T when the pulse 
has an FWHM duration t « 0.42 T. 
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